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The  aim  of  this  project  initially  was  to  develop  accurate  non-parametric 
methods  for  determining  the  energy  band  structure  of  semi -conducting  polymers 
and  then  to  incorporate  within  the  same  model  a  method  to  determine  the  elec¬ 
tron  levels  due  to  impurity  atoms  used  as  replacements  for  the  side  groups 
about  the  polymer  back  bone.  Secondly,  we  were  to  study  the  energy  structure 
of  several  polymers  with  and  without  impurity  side  groups.  These  tasks  have 
been  successfully  accomplished . 

The  essential  approximation  employed  here  is  that  of  Hartree-Fock.  Here 
the  many  body  Hamiltonian  is  simply: 


-i  K- 1 

i=l  i,I 


I 


+  7  l 


i,j  i ri-f j 


+  I  f  — O— 
2  IMjl 


The  trial  wavef unction  has  the  form: 


?(?!... rHJ  -  7=y  detj4.iCr1)  |  .  (2) 

O 

Here  we  use  atomic  units,  e=1i  =  m=l,  length  in  a.u.  (1  a.u.  =  .53A)  and 
energy  in  Hartree's  (1  Hy=27.2  eV)  .  Lower  case  letters  refer  to  electron 
properties  and  upper  case  to  nuclear  properties.  If  we  constrain  the  <f>'s  to 
be  orthonormal  and  eigenstates  of  s z  and  require  that  the  expectation  value 
of  H  be  minimized  with  respect  to  our  trial  ty,  then  the  <f>'s  are  determined 
by  the  unrestricted  Hartree-Fock  eq.: 


P(rr')  . 


(4) 
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#(r r')  is  the  operator  which  interchanges  coordinates  r  and  r ' .  It  is  this 

set  of  equations  which  we  will  solve.  The  final  results  will  have  one  well 

characterized  deficiency.  The  usual  Hartree-Fock  solution  has  too  great  a 

gap  between  the  uppermost  occupied  orbital  (HOMO)  and  the  least  unoccupied 

virtual  orbital  (LUMO) .  This  deficiency  has  been  well  studied  in  both  polymer 

and  solid  state  systems^  and  is  only  of  concern  when  one  needs  the  value  of 

the  gap  precisely.  In  any  event  this  deficiency  may  be  readily  corrected  by 

2/ 

using  the  electronic  polaron  for  long  range  correlation—  and  AE(SCF)  tech- 

3/ 

niques  for  short  range  orbital  relaxation—  .  The  internal  band  details  are 
far  less  sensitive  to  such  corrections.  For  the  express  purposes  of  this 
report  we  ignore  such  corrections  other  than  to  say  their  inclusion  is  quite 
trivial . 

The  essential  problems  in  solving  the  system  of  equation  (4)  lie  in  that 
some  simplifying  assumptions  are  needed  still,  and  it  is  with  these  that  we 
initially  concern  ourselves.  To  see  the  problems  and  procedure,  consider  as 
a  prototype  system  polyethylene,  which  has  a  unit: 

H  H 

I  _  I 

—  C  C 

H  H 

This  unit  is  repeated  infinitely.  If  one  uses  energy  band  techniques,  this 
repetition  must  have  one  dimensional  translational  symmetry  and  hence  one 
must  idealize  out  of  existence  the  random  twisting  of  the  backbone  of  such  a 
polymer  which  occurs  in  the  real  world.  Additionally,  in  the  infinite  system 
Ewald  sums  occuring  in  F  yield  a  poorly  defined  energy  zero  and  thus  accurate 
positioning  of  Impurity  levels  without  use  of  methods  such  as  Stater-Koster 


becomes  hard . 


There  is  an  alternate  approach.  This  is  to  reduce  the  polymer  to  a 
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finite  size.  This  is  the  cluster  problem.—  Here  the  essential  approxima¬ 
tion  is  to  reduce  an  infinite  size  system  to  a  small,  finite  one.  In  doing 
this  one  loses  the  e(k)  curves,  beloved  of  the  band  theorist,  but  retains  a 
set  of  energy  levels  which  have  no  k  label.  These  levels  form  a  state  density 
which  is  consistent  with  those  of  energy  band  theory  provided  one:  a)  uses 
proper  boundary  conditions  on  the  cluster^;  b)  uses  a  large  enough  cluster.—^ 
Tue  use  of  clusters  compensates  for  these  problems  by  permitting  one  to  study 
non-ideal  geometries,  and  by  allowing  simple  inclusion  of  defects  or  impuri¬ 
ties  in  the  calculation.  We  tried  both  methods  initially  and  decided  finally 
that  the  cluster  techniques  were  to  be  preferred  for  this  study.  In  the 
remainder  of  this  report  we  concentrate  on  cluster  results. 

The  size  of  cluster  has  been  investigated  by  J.-M.  Andre  and  J.  Delhalle 
for  the  case  of  polyethylene  using  the  Extended  Huckel  method.—^  They  study 
the  energy  levels  for  clusters  of  1,  2,  3,  4,  5  and  6  C-atoms  and  for  the 
idealized  band  geometry  ( ® C-atoms) .  Essentially,  the  results  are  fully 
converged  in  band  position  with  1  C-atom  and  in  width  with  4  C-atoms.  In  our 
study  we  use  6  C-atoms  in  the  backbone  as  a  guarantee  of  convergence.  The 
second  problem  is  the  finite  cluster  termination.  Polyethylene  is  a  covalent 
system  and  according  to  the  prescription  of  Kunz  and  Klein  one  should  place  an 
excess  electron  and  nucleus  at  each  end  of  the  intercepted  chain.— ^  Conven- 

4/ 

tionally  one  chooses  a  proton  for  the  replacement  nucleon—  ,  but  it  has  been 
observed  by  Allen  that  one  may  vary  the  charge  of  this  nucleon  to  force  the 
atoms  at  the  end  of  the  chain  to  have  core  properties  similar  to  those  in  the 
center.—^  We  concur  with  this  suggestion  and  used  it  here.  We  found,  however, 
when  the  charge  was  varied  in  units  of  0.1  electron  charges  for  our  system  that 
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a  value  of  1  (the  proton  charge)  was  in  fact  optimal.  We  believe  this  Is 
consistent  with  polyethylene  being  an  essentially  pure  covalent  system. 

The  cluster  energy  levels  for  polyethylene  are  shown  in  Figure  1  for 

the  valence  levels  as  given  in  our  calculation.  In  addition  the  density  of 

States  for  polyethylene  in  all  available  alternate  models  using  band  theory 

are  shown  in  Figure  1  as  well  by  means  of  comparison.  It  is  clear  that 

a  reasonable  agreement  exists  between  our  UHF  results  and  the  existing  ab 
8  / 

initio  result.—  Agreement  with  the  other  results  all  of  which  are  either 

9/ 

highly  approximate  or  parametric  is  less  good.  These  include  F.S.G.O.—  , 
L.C.L.O.— f ,  E.H.C.O.— ,  CNDO/2^  and  CNDO/2^  methods.  The  disagreement, 
small  though  it  is,  between  our  UHF  results  and  the  ab  initio  results  probably 
reflects  the  level  of  Hartree-Fock  achieved  in  the  ab  initio  results  and  the 
integral  approximations  used  rather  than  the  finite  cluster  size.  We  use  a 
gaussian  basis  set  of  double  zeta  accuracy  here  and  evaluate  all  our  integrals 
exactly  to  the  precision  permitted  by  our  computer  (8.5  significant  figures 
here) .  We  conclude  that  our  technique  produces  a  viable  energy  level  scheme 
for  the  polymer  of  interest  and  is  competitive  in  this  respect  to  energy  band 
calculations.  A  like  study  was  performed  for  a  SI  based  polymer  in  a  hypothe¬ 
tical  structure  similar  to  polyethylene.  Results  were  quite  similar. 

The  final  study  concerned  the  doping  of  the  polymer  chain  with  impurity. 

In  this  event  a  hydrogen  was  removed  from  an  internal  C(Si)  and  an  impurity 
substituted.  In  the  case  of  [CHj]^  wt  used  Li,  Be,  B  and  F  and  for  [SiH]^ 
we  used  Li,  Be  and  F.  We  attempted  to  see  if  for  polyvalent  impurities,  in 
this  case  Be  on  [ScH]„,  one  could  remove  a  second  H  atom  as  well.  We  find 
this  to  Indeed  be  a  lower  energy  state  however,  the  one  electron  energy  levels 
are  not  much  modified  by  this  and  we  ignore  this  complexity  in  the  other  studies 
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The  one  electron  H.F.  levels  for  [CH^]^  and  with  Li,  Be,  B,  and  F  are  seen 
in  figure  2  and  for  [SiH^ljj  and  with  Li,  Be  and  F  in  figure  3.  It  is  clear 
from  these  figures  that  the  effect  of  the  impurity  is  to  allow  the  occupied 
levels  to  remain  largely  undisturbed,  except  that  occupied  levels  are  now 
found  in  the  energy  gap.  The  effect  on  the  virtual  levels  is  far  more  pro¬ 
nounced  and  causes  the  formation  of  local  trap  states  below  the  conduction 
band.  The  apparent  number  of  these  is  much  exaggerated  here  due  to  the  finite 
size  of  the  cluster.  In  actuality,  the  electronic  structure  would  look  like 
a  superposition  of  the  perfect  polymer  structure  and  that  shown  including 
impurity.  Finally,  we  observe  here  that  including  correlation  effects  would 
reduce  the  band  gap  by  about  0.2  Hy.  Thus,  we  further  conclude  that  proper 
modification  of  side  groups  in  the  semiconducting  polymers  offers  a  viable 
means  of  gap  modification  but  also  introduces  the  hazard  of  trapping  states 
in  the  gap.  In  this  respect  such  systems  behave  much  as  ordinary  semiconductors. 

The  results  obtained  are  being  combined  into  form  to  submit  for  publication. 
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Figure  1:  The  one  electron  occupied  valence  levels  of  polyethylene  from  our 
study  are  shown  in  a.  In  b  we  find  a  state  density  for  this 
system  from  ab  initio  energy  band  theory,  in  c_  we  see  like  results 
using  a  FSGO  method,  in  d^  using  a  LCLO  method,  in  e_  using  the  EHCO 
method,  in  f_ using  the  CNDO/2  method  and  in  using  a  CNDO/w  method. 


Figure  2:  The  one  electron  energy  levels,  both  occupied  and  virtual  for  poly¬ 
ethylene,  including  impurities  are  seen  for  several  cases.  The 
fermi  energy,  e^,  is  indicated. 


Figure  3:  The  one  electron  energy  levels,  both  occupied  and  virtual  for  [Sil^] , 
including  impurities  are  seen  for  several  cases.  The  fermi  energy, 
is  shown. 
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